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1. Introduction 



The infrared behavior of Green's functions in Landau gauge has been the topic of numerous 
lattice studies by several groups in the past few years. Particular attention has been devoted to 
the gluon and ghost propagators, whose infrared behavior is at the heart of the Gribov-Zwanziger 
confinement scenario [0, ]3L Q]. There is now a consistent picture — from extensive numerical 
simulations on very large lattices [§, 0, 10, 11, [l2|, 13, [l4]] — that (in three and in four space- 
time dimensions) the Landau gluon propagator shows a massive solution at small momenta and 
that the Landau ghost propagator is essentially free in the same limit. These results are not in 
agreement with the original Gribov-Zwanziger scenario [[|, [|] but they can be explained in the so- 
called refined Gribov-Zwanziger framework [O]. Let us also recall that a massive gluon allows a 
better description of experimental data Jl6| ] and it has been related to color confinement by various 



authors [17, 18]. 



Since the evaluation of Green's functions depends on the gauge condition, it is important 
to consider different gauges in order to obtain a clear (possibly gauge-independent) picture of 
color confinement. Needless to say, this investigation should be carried out at the nonperturbative 
level. This is done from first principles using lattice simulations. In addition to the Landau gauge 
case cited above, numerical studies of Green's functions have also been done in Coulomb gauge 



J19| , pQ| , pT| ], A -gauge (a gauge that interpolates between Landau and Coulomb) y22| ] and maximally 



Abelian gauge [J23J, [24J]. For interesting recent comparisons of results in Landau and in Coulomb 
gauge see [|25], [2q ]. 

On the other hand, the linear covariant gauge — which is a generalization of Landau gauge — 



proved for a long time quite hostile to the lattice approach [27, 28, 29, 3(| 31, 32, 33|, [34|]. Recently 
we have introduced a new implementation of the linear covariant gauge on the lattice |JXJJ , based 
on a minimizing functional that extends in a natural way the Landau case while preserving all the 
properties of the continuum formulation. Let us note that, having a minimizing functional for the 
linear covariant gauge allows a numerical investigation of the first Gribov region Q. for the case 
of gauge parameter B, ^ 0. Such an investigation has been done analytically in for a small 
value of Bj , but a similar numerical study is still lacking. At the same, a numerical investigation 
of the infrared behavior of gluon and ghost propagators at B, / could provide important inputs 
for analytic studies based on Dyson-Schwinger equations [|(], Finally, it has been recently 
that the background-field Feynman gauge is equivalent (to all orders) to the 



proven [Bq, 39 



pinch technique [f40|, |4l|]. Thus, numerical studies using the Feynman gauge, which corresponds 
to the value B, = 1, will allow a nonperturbative evaluation of the gauge-invariant off-shell Green's 



functions of the pinch technique [42]. 



2. Linear Covariant Gauge 



In the linear covariant gauge the gluon field A* (x) satisfies (in the continuum) the relation 



d^{x) = A h (x), 



(2.1) 
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where A b (x) are real- valued functions generated using a Gaussian distribution 

1 



A b (x) 



expl-^^[A b (x) 



(2.2) 



with width yf~^. 

The limit E, — > corresponds to the standard Landau gauge. In this case, the gauge condi- 
tion is (classically) equivalent to the Lorenz-gauge (sometimes mistakenly called Lorentz-gauge) 
condition [p3l 

. (2.3) 



d„A b Jx) 



This condition can be imposed by minimizing the functional 

<?w{A g } - / d 4 x^\(A%(x) 



(2.4) 



with respect to the gauge transformations {g(jc)}. Let us recall here that, from the second variation 
of (^lgIA 8 }, we can define the Faddeev-Popov operator Then, for the gauge-fixed configura- 
tions, i.e. for local minima of <^g{A s }, we have that this operator is positive-definite. This set of 
local minima defines the first Gribov region Q. ^ [|] . 



In Ref. [ p7[ ] it was shown that a similar minimizing functional S'lcgIA 8 } for the linear covari- 
ant gauge — i.e. for § / — does not exist. Indeed, if it existed, we could write 



<%lcg [A 8 , A] = ^g[A«] + &\A*,A] , 



(2.5) 



for some functional J£"[A S ,A]. Then, the second variation of S'lcg with respect to the gauge trans- 
formation g(x) = e w W would satisfy the relation 



d 2 & 



LCG 



•LCG 



dw b {x)dw c {y) dw c (y)dw b {x) 



(2.6) 



On the other hand, one can show that these two terms are, respectively, proportional to the structure 
functions f acb and f abc . Since these functions are completely anti-symmetric in the color indices, 



this equality cannot be realized [27] 



3. Bypassing the No-Go Condition 

One can of course avoid the no-go condition above by relaxing its hypotheses. For example, 
one can consider a different gauge condition, such as 



F[d,A b (x)-A b (x)] = 



(3.1) 



with F[0] = 0, for which a minimizing functional exists. Indeed, it has been shown in Ref. pJ\\ 
that the minimizing functional 



■' Ll.b I 



d 4 x£\ d^W-A'ix) 



(3.2) 
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allows one to impose the gauge condition 



where D" h is the covariant derivative. 



■A b (x) 



0, 



(3.3) 



On the other hand, the use of a different gauge condition introduces new problems [27]. For 
example, one can bring in spurious solutions, corresponding to F[s] = for s ^ 0. In the above 
case, these solutions are the zeros of the operator D" b d v . Also, in general, the second variation 
of the minimizing functional, or equivalently the first variation of the gauge-fixing condition (fO|), 
does not correspond to the Faddeev-Popov operator ^# = —d^D ab of the usual linear covariant 
gauge. Finally, the lattice discretization of the functional (3^2) is not linear in the gauge transfor- 
mation {g(x)}. This makes the numerical minimization difficult and one has to rely on a specific 
discretization of the minimizing functional [28, ^] in order to make the lattice approach feasible. 

More recently [ |33| ] the no-go condition has been overcome by avoiding the use of a minimizing 
functional <§wG{A g }- To this end, following the perturbative definition of the linear covariant gauge 
in the continuum, one first fixes the gluon field to Landau gauge d^A 1 ^ (x) = 0. Then, one considers 
the equation 

(3.4) 



d^f) (x) = A b (x) . 

The solution <p c (x) of this equation can be used as a generator of a second gauge transformation. 
Note that, after fixing the lattice (or minimal) Landau gauge, the operator —d^D^' is positive- 
definite and can be easily inverted. For small <p c (x), the final (gauge-transformed) gluon field 
A' (jc) satisfies the condition 



d^A'fa) 



^K+D^')(x) 



A b (x) 



(3.5) 



Of course, the above result is correct only for infinitesimal gauge transformations. On the other 



hand, usually (f) c (x) is not small in a numerical simulation. Indeed, numerical tests [33] have shown 



that the distributions of d^A' (x) and of A (x) do not agree very well. At the same time, the relation 
p 2 Di(p 2 ) = valid in the linear covariant gauge for the longitudinal gluon propagator D[(p 2 ), is 



also not well verified by the data at small momenta [33]. 



4. A New Approach 

Our new approach [|I|] is based on removing an implicit hypothesis of the no-go condition, 
i.e. that the gauge transformation {#(*)} appears in the minimizing functional in the "canonical" 
way A g . Thus, we may look for a minimizing functional of the type <%wg{A- s ,g} instead of simply 
(?lcg{A 8 }- Indeed, the lattice linear covariant gauge condition can be obtained by minimizing the 
functional 1 

<?Lcc{U g ,g} = £lg{U 8 } + 91 7>£ ig(x)A(x) , (4.1) 

x 

where 

<?L G {U 8 } = -XTrY, g(x)U fl (x)g^x + e^ (4.2) 
x,n 

'Note that solving a system of equations By = £ is equivalent to minimizing the quadratic form jV/Si// — \jf£. 
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Figure 1: Convergence of the numerical gauge fixing. Here we report the value of A = S^j,[V • A b {x) — 
A h (x)] 2 as a function of the number of iterations n for a given configuration. Left: j3 = 4, V = 8 4 , E, = (red 
line) and E, — 0.5 (green line). Right: j3 = 4, V = 16 4 , £ = (red line) and <^ = 0.05 ( green line). Note the 
logarithmic scale on the y axis. 



is the minimizing functional for the lattice Landau gauge. Here, the link variables (x) and the 
site variable g(x) are matrices belonging to the S\J(N C ) group (in the fundamental representation). 
We also indicate with 9? the real part of a complex number and with Tr the trace in color space. 
Note that the functional $lcg{U s ,g} is linear in the gauge transformation {g(x)}. 

By considering a one-parameter subgroup g(x, t) = exp [/Ty*(x)A*] of the gauge transforma- 
tion {g(x)} it is easy to check that the stationarity condition 

= V y*(jc) (4.3) 

implies the lattice linear covariant gauge condition 

V-A b (x) = £A£(x) - Aj(*-e M ) = A b (x) . (4.4) 

Here, X b are the traceless Hermitian generators of the Lie algebra of the SU(A^ C ) gauge group, 
satisfying the usual normalization condition 

7Y(VA C ) =28 bc . (4.5) 

Also, we used A b (x) = Tr[A(x) X b \ and, similarly, Aj(jc) = Tr[A^(x)X b ] . 

At the same time, the second variation (with respect to the parameter %) of the term ig(x) A(x) , 
on the r.h.s. of Eq. (fO|), is purely imaginary. Thus, it does not contribute to the second variation of 
the functional <%lcg{U 8 ,g}- This implies that, using the above minimizing functional, one finds for 
the Faddeev-Popov matrix ^# a discretized version of the usual Faddeev-Popov operator —d ■ D. 



5 



Lattice linear covariant gauge 



Attilio Cucchieri 



5. Numerical Tests 



We have performed some numerical tests 2 with the functional (11), using the stochastic- 
overrelaxation algorithm [[P], ^5, 46]. 



For these first tests we considered the 4d SU(2) case at /3 = 4, for V = 8 4 and 16 4 , with 
% = 0.01,0.05,0.1 and 0.5. The numerical gauge fixing works very well when E, ^ 0, at least for 
relatively small lattice volumes V and gauge parameter t, (see the next section). In Figure [l] (left 
panel) we compare the gauge fixing for a given configuration at j3 = 4 and V = 8 4 for the Landau 
case % = (red line) and for £ = 0.5 (green line). The rate of convergence is essentially the same in 
the two cases. However, note that the tuning of the stochastic-overrelaxation algorithm is different. 
Indeed, we found that, for the chosen configuration, the best value for the stochastic-overrelaxation 
parameter p was about 0.73 in the Landau case and about 0.5 when £ = 0.5. A similar result (see 
Figure [T], right panel) is obtained for j8 = 4 with V = 16 4 for the Landau case (E, = 0, red line) 
and for £ = 0.05 (green line). In this case the best value for p was about 0.82 in the Landau case 
and about 0.81 when % = 0.05. Let us note that the functional (§ > LCc{U g ,g} can be interpreted as 
a spin-glass Hamiltonian for the spin variables g(x) with a random interaction given by U^(x), in 
a random external magnetic field A(x). The presence of this magnetic field does not modify the 



convergence matrix [47] or, as a consequence, the behavior of the algorithm. 3 

We also checked that the quantity p 2 Di{p 2 ) is constant within statistical fluctuations in all 
cases considered. For V = 16 4 and t, = 0. 1 and 0.5 the data are shown in Figure |2| In these cases, 
a fit of the type a/p b for Di(p 2 ) gives a = 0.0994(7), b = 2.003(9) with a % 2 /dof = 0.9 when 
% =0.1 and a = 0.502(5), b = 2.01(1) with a % 2 /dof = 1.1 when % =0.5. Similar fits have been 
obtained in the other cases. 



6. Discretization of the Gluon Field 



In the above tests we used the usual discretization 



A M (x) 



■Ufa) 



traceless 



2i 



(6.1) 



for the gluon field. However, one has to recall that, using this standard (compact) discretization, 



the gluon field is bounded. On the other hand, the functions A b (x) [see Eqs. (2T ) and (2.2)] satisfy 
a Gaussian distribution, i.e. they are unbounded. This can give rise to convergence problems [|^]. 
Moreover, the problem is more severe for a larger width of the Gaussian distribution. 

A possible way out of this problem is, of course, the use of different discretizations of the 
gluon field A^(x), in order to improve the convergence of the minimizing algorithms. To this end 
we also did some tests using the angle projection [50] and the stereographic projection [|5l|]. Note 
that, in the last case, the gluon field is in principle unbounded even for a finite lattice spacing. We 



2 Note that Eq. ( |4.4| ) implies £ t A b (x) = 0. Thus, after generating the functions A b (x) using the Gaussian distribution 
(2.2), one has to remove possible zero modes from them. 



3 Note that, in the case j8 = °° p^], the Landau case corresponds to solving the Laplace equation while the linear 
covariant gauge corresponds to solving a Poisson equation. As is well-known, these two equations show the same critical 
behavior when solved using a relaxation method such as Gauss-Seidel mm. 
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Q 0.52 




Figure 2: The longitudinal gluon dressing function p 2 D/(p 2 ) as a function of the lattice momentum p (in 
lattice units). We also show the predicted value p 2 Di(p 2 ) = %. Left: j3 = 4, V = 16 4 and % =0.1. Right: 
/3 = 4, V = 16 4 and t, =0.5. Note the relatively small range of values on the y axis. 



Table 1: Smallest value of j3 for which the numerical gauge-fixing algorithm showed convergence for the 
lattice volume V — 8 4 . Results are reported for the three different discretizations considered and for five 
different values of the gauge parameter % . For each case we used five different configurations. 



% 


stand, disc. 


angle proj. 


stereog. proj. 


0.01 


2.2 


2.2 


2.2 


0.05 


2.2 


2.2 


2.2 


0.1 


2.2 


2.2 


2.2 


0.5 


2.8 


2.6 


2.5 


1.0 




3.0 


2.5 



also stress that, in both cases, the numerical implementation gets simplified if one uses the Cornell 
method 46] instead of the stochastic-overrelaxation algorithm, as done in the simulations 

reported in the previous section. 

We tested the standard discretization, the angle projection and the stereographic projection 
using V = 8 4 , § =0.01,0.05,0.1,0.5,1.0 and J3 = 2.2,2.3, ... ,2.9,3.0. We found (see Table |) 
that the stereographic projection allows one to simulate at slightly larger values of t, , for a given 
lattice volume V and lattice coupling j3, compared to the other two cases. 4 



7. Continuum Limit 

In Ref. [ p3[ ] it was shown that in the SU(A^ C ) case, in order to obtain the correct continuum limit, 
the functions A b (x) should be generated using a Gaussian distribution with width y/a = ^2N C £, /j8, 

4 Note that, as explained in the next section, the Gaussian distribution generated in the simulation actually has 
squared width cr = 4|/j3 [in the SU(2) case]. 
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instead of the width y£ . Thus, for /3 < 2N C the lattice width yfa is even larger than the continuum 
width On the other hand, one can always obtain a sufficiently small value for a by considering 
large enough values of the lattice coupling /?. However, if /3 is too large, the physical volume is too 
small (for a given lattice size) and one cannot study the infrared limit of the theory. 

Note that, in the SU(2) case, one has a = £, only for j8 = 4, which corresponds to a very small 
lattice spacing, i.e. a w 0.001 fm [52]. On the contrary, in the SU3 case, one has a = % for j8 = 6, 
corresponding to a lattice spacing a = 0.102 fm [53], usually employed in lattice numerical studies. 
Thus, simulations for the linear covariant gauge are probably easier in the SU(3) case. 

One should also recall that the gluon field A and the gauge parameter % are (multiplicatively) 
renormalized by the same factor Z3, i.e. Ab = Z 3 Ar and £,b = Z3 Qj. Here, B and R indicate bare 
and renormalized quantities respectively. This implies that, on the lattice, data obtained for two 
different values of j8 in the scaling region — e.g. j3i and p 2 — would give the same (renormalized) 
propagator only if the multiplicative factor Rz = (j3i ) / (J82 ) relating the propagators 5 also re- 
lates the gauge parameters %\ and <§2- Since the value of Rz is not known a priori, one has to find 
numerically pairs of parameters Q3 , £ ) yielding the same continuum renormalized propagators. 



8. Transverse Gluon Propagator 

In this section we present preliminary results for the momentum-space transverse gluon prop- 
agator D t (p 2 ) for different values of Using the stereographic projection, we have simulated 
the linear covariant gauge at j3 = 2.2 and j8 = 2.3 for the lattice volumes V = 8 4 , 16 4 and 24 4 , 
considering several values of the gauge parameter E, in the SU(2) case (see Table |2|). We observe 
that the quantity Di(p 2 )p 2 /£,, which should be equal to 1, has a value of 0.999(2) when averaged 
over all data D/(p 2 ) produced. From the results shown in Figures [3] and ^ one clearly sees that, 
as in Landau gauge, the propagator is more infrared suppressed when the lattice volume increases. 
At the same time, for a fixed volume V, the propagator is also more infrared suppressed when the 



gauge parameter % increases. The latter result is in agreement with Ref. jplQ. 

Let us recall that, on the lattice, the finite size of the system corresponds to an infrared cutoff 
~ 2n/L, where L is the lattice size. In the four-dimensional case for j8 = 2.2, in Landau gauge, one 
needs to use a lattice volume V of about 64 4 in order to obtain infinite- volume-limit results [Q, [It]]. 
The data shown in Figures || and |] seem to indicate that similar lattice volumes are also needed 
for the linear covariant gauge. On the other hand, considering the data reported in Section ^, the 
extrapolation to infinite volume for a given /3 and a fixed value of E, seems harder than in Landau 
gauge. Indeed, as V — > °°, the number of sites characterized by a large value for the function A b (x) 
increases, making the convergence of the gauge-fixing method more difficult. 



9. Conclusions 

We have discussed a new lattice implementation of the linear covariant gauge. The gauge 
fixing is done using a minimizing functional S'lcciU 8 ,g}, which is a simple generalization of the 
Landau-gauge functional $lg{U 8 }- Tests done for the SU(2) case in four space-time dimensions 

5 See Ref. p4[] for the case of Landau gauge. 
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Table 2: Values of the lattice coupling /3, the lattice volume V and the gauge parameter S, used in our 
simulations. We also report, in each case, the total number of configurations considered. 



p 


V 




number of conf. 


j8 


V 


% 


number of conf. 


2.2 


8 4 


0.0 


500 


2.3 


8 4 


0.0 


500 


2.2 


8 4 


O.Ol 


500 


2.3 


8 4 


0.01 


500 


2.2 


8 4 


0.05 


500 


2.3 


8 4 


0.05 


500 


2.2 


8 4 


O.l 


500 


2.3 


8 4 


0.1 


500 


2.2 


8 4 


0.2 


500 


2.3 


8 4 


0.3 


500 


2.2 


8 4 


0.3 


500 


2.3 


8 4 


0.5 


500 


2.2 


_ A 

8 4 


0.4 


500 


2.3 


8 4 


0.6 


500 


2.2 


_ A 

8 4 


0.5 


542 


2.3 


8 4 


0.7 


543 


2.2 


16 4 


0.0 


400 


2.3 


16 4 


0.0 


400 


2.2 


16 4 


0.01 


400 


2.3 


16 4 


0.01 


400 


2.2 


16 4 


0.05 


400 


2.3 


16 4 


0.05 


400 


2.2 


16 4 


0.1 


400 


2.3 


16 4 


0.1 


400 










2.3 


16 4 


0.2 


336 


2.2 


24 4 


0.0 


158 


2.3 


24 4 


0.0 


200 


2.2 


24 4 


0.01 


200 


2.3 


24 4 


0.01 


200 


2.2 


24 4 


0.05 


160 


2.3 


24 4 


0.05 


310 










2.3 


24 4 


0.07 


59 



show that this approach solves most problems encountered in earlier implementations and ensures 
a good quality for the gauge fixing with a ratio Di(p 2 )p 2 « 1 for all cases considered. We have 
also presented preliminary results for the transverse gluon propagator D t (p ). 

As discussed in Sections || and |g, the only open problem is the convergence of the gauge- 
fixing algorithm at large lattice volumes when the gauge parameter £ is also large. However, as 
mentioned above, this problem is probably less severe for the SU(3) group compared to the SU(2) 
case. We are currently simulating other values of /3 and E, in the 4d SU(2) case and considering 
simulations also of the SU(3) group and of the 3d case J55| ] . 
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